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On the basis of a next-to-leading-order calculation in chiral perturbation theory, the first complete 
analysis of isospin breaking for direct CP violation in K° — > 2n decays is performed. We find a 
destructive interference between three different sources of isospin violation in the CP violation 
parameter e'. Within the uncertainties of large- N c estimates for the low-energy constants, the 
isospin violating correction for e is below 15 %. 



Direct CP violation in K° — > 2tt decays is now estab- 
lished to an accuracy of less than 10 %: 



Re e'/e = (1.66 ±0.16) x 10" 



(1) 



The theoretical situation is less satisfactory due to large 
uncertainties in the hadronic matrix elements of four- 
quark operators. Although several calculations exist that 
come close to the experimental result (e.g., Ref. 0) the 
accuracy is still rather limited. 

K decays involve a delicate interplay between electro- 
weak and strong interactions in the confinement regime. 
Chiral perturbation theory (CHPT) provides a conve- 
nient framework for a systematic low-energy expansion 
of the relevant amplitudes. In this letter, we concentrate 
on a quantitative analysis of isospin violation in the CP- 
violating parameter e'. The results presented here are 
part of a complete calculation of K — > 27r decays to 
next-to-leading order in the chiral expansion and to first 
order in isospin violation, including both strong isospin 
violation (m u ^ m^) and electromagnetic corrections. 

Although isospin violation is in general a small effect 
it must be included 0, 0, IE m a precision calculation 
of e' because it affects the destructive interference be- 
tween the two main contributions to e' from normal and 
electromagnetic penguin operators. 

The main features of our approach are the following: 

1. We include for the first time both strong and elec- 
tromagnetic isospin violation in a joint analysis. 

2. Nonleptonic weak amplitudes in CHPT depend on 
a number of low-energy constants (LECs) : we use leading 
large- N c estimates for those constants 0, El • Uncertain- 
ties arise from (i) input parameters in the leading 1/N C 
expressions as well as from (ii) potentially large sublead- 
ing effects in l/N c . We discuss the impact of both (i) 
and (ii) on the relevant quantities. 

We adopt the following parametrization of K — > nir 
amplitudes 
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In the limit of CP conservation, the amplitudes Aq,A 2 , 
and A~2 are real and positive. In the isospin limit, A 2 — 
A 2 j X2 = xt m the standard model and the phases \i 
coincide with the corresponding pion-pion phase shifts at 
-Kerns = Mr- In terms of these amplitudes, the direct CP 
violation parameter e' is given by 
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For the systematic low-energy expansion of K — > ttt: 
amplitudes to next-to-leading order, a number of effective 
chiral Lagrangians is needed. Here, we only write down 
the relevant nonleptonic weak Lagrangian, referring to 
Ref. |sj for further details: 



Avcak = G 8 F 4 (L^) 23 + G 8 g cwk e 2 F 6 (U^QU) 
+ G 27 F A (l M 23L^ + \ 
+ GsN^Of + Y, G 27 D t F 2 Of + h.c 



: -£^2 l^i3 
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The matrix L M = iWD^U represents the octet of V — A 
currents to lowest order in derivatives where the SU(3) 
matrix field U contains the pseudoscalar fields. Q rep- 
resents the light quark charge matrix. In addition to 
the lowest-order couplings G%,G 2 -j, and Gs^ewk, this 
Lagrangian contains LECs G%Ni, G 27 Di of 0(GfP 4 ) 
0. Due to the disparity in size between G$ and G 27 
(AI =1/2 rule), we consider isospin violating effects 
in the octet part only. More precisely, we work up 
to 0{G s p\G 8 (m u - m d )p 2 ,G 8 e 2 p 2 ) and to 0{G 27 p 4 ) 
for octet and 27-plet amplitudes, respectively. For this 
purpose, the chiral effective Lagrangian for strong in- 
teractions to 0(p 4 ), the electromagnetic Lagrangian to 
0(e 2 p 2 ) and the electroweak Lagrangian to 0(G%e 2 p 2 ) 
are also needed 13J. 
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The resulting decay amplitudes depend on a large num- 
ber of renormalized LECs. In addition to the well-known 
strong LECs of 0(p 4 ), we make use of existing estimates 
for the electromagnetic couplings of 0(e 2 p 2 ) [9j. For the 
nonleptonic weak LECs of 0(GfP 4 ) and 0(G$e 2 p 2 ), we 
employ leading large- N c estimates 0, The starting 
point is the effective AS = 1 Hamiltonian in the three- 
flavour theory, 

( '" " "* ^a(/iSD)Qi(/iSD) , (5) 
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H eff 



-V jV* , . 

^2 ud us / j 



obtained after integrating out all fields with masses larger 
than /xsD — 1 GeV. The local four-quark operators Qi 
contain only the light degrees of freedom whereas the 
Wilson coefficients C; are functions of heavy masses and 
CKM parameters [3 EH: 

Ci (MSd) = Zi(fisu) + TJ/i(/i SD ) 



V ., V* 
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All CP-violating quantities are proportional to Im r. In 
our analysis, only ratios of CP-violating amplitudes will 
appear so we do not need an explicit value for r. 

At leading order in l/N c , the matching between the 
three-flavour quark theory and CHPT can be done ex- 
actly because the T-product of two colour-singlet quark 
currents factorizes. Since quark currents have well-known 
realizations in CHPT, the hadronization of the weak op- 
erators Qi is straightforward [2J, y|. The weak LECs 
are expressed in terms of strong LECs of order p 4 , p 6 , 
and e 2 p 2 . The 0(p 4 ) and 0(p 6 ) couplings arising from 
bosonization of Qg are estimated within resonance sat- 
uration and are determined by scalar exchange (for de- 
tails see Refs. |3E3)- 

Although admittedly this is not 
a complete estimate for the 0(p 6 ) couplings it certainly 
provides the correct order of magnitude. 

Subleading effects in 1/N C are known to be sizable at 
leading chiral order. We will therefore not use the large- 
N c values for ReGs, ReG 2 7 in the numerical analysis 
but instead determine these couplings from a fit to the 
K — > 7T7T branching ratios 31- The other combination 
of interest is the ratio Im(G8<7ewk)/ImG8. In this case 
the size of 1/N C effects is approximately determined by 
existing calculations going beyond factorization [13|. It 
turns out that corrections to this ratio are not very large 
(roughly —30%). The dominant uncertainty comes then 
from the input parameters in the factorized expressions. 

At next-to-leading chiral order we estimate subleading 
effects in 1/N C by varying the chiral renormalization scale 
at which the large- N c results are supposed to apply. We 
consider two options for the higher-order couplings G$Ni, 
. . . in the Lagrangian jljl: one may either adopt the large- 
N c predictions for the ratios (GgN^/Gg or directly for 
the couplings GsNi. The uncertainties related to these 
two choices will be taken into account. We also include 
in the errors the effect of changing the short-distance 
renormalization scale /xsd between M p and 1.3 GeV, with 
the central value at 1 GeV. 



Whenever the errors are sizable, they are however 
largely dominated by the intrinsic uncertainty (an ef- 
fect of higher order in 1/N C ) at which chiral scale 
fCH the large-TV^ estimates for the renormalized LECs 
GgN( (ugh)) ■ ■ ■ apply- We vary this scale between 0.6 
GeV and 1 GeV, with the central values at M p . 

We now turn to the different sources of isospin viola- 
tion in e'. The expression © is valid to first order in CP 
violation. Since ImA/ is CP odd the quantities ReA/ and 
Xi are only needed in the CP limit (/ = 0, 2). We disre- 
gard the phase that can be obtained from the K — > irir 
branching ratios. The same branching ratios are usually 
employed to extract the ratio ujs = ReA 2 /ReAo assum- 
ing isospin conservation. Accounting for isospin violation 
via the general parametrization J5J), one is then really 
evaluating lo + = ReA^ /ReAo rather than wg. The two 
differ by a pure AI = 5/2 effect: 
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Because o;+ is directly related to branching ratios it 
proves useful to keep o>+ in the normalization of e', in- 
troducing the AI = 5/2 correction / 5 / 2 

Since Irm42 is already first order in isospin violation 
the formula for e' takes the following form, with all first- 
order isospin breaking corrections made explicit: 
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where 
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and the superscript (0) denotes the isospin limit. 

To isolate the isospin breaking corrections in e', we 
write the amplitudes Aq , A 2 more explicitly as 
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where <5Ai/ 2j 3/ 2 , "4-5/2 are fast order in isospin violation. 
To the order we are considering, the amplitudes Aai have 
both absorptive and dispersive parts. To disentangle the 
(CP conserving) phases generated by the loop amplitudes 
from the CP-violating phases of the various LECs, we 
express our results explicitly in terms of T>isp Aai and 
Abs Aai- 

To first order both in CP violation and in isospin break- 
ing, we obtain 
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where \Ai?\ = ^ {Re[Visp A^]) 2 + (Re As A^]) 2 . 

These expressions are general results to first order in 
CP and isospin violation but they are independent of the 
chiral expansion. Working strictly to a specific chiral or- 
der, these formulas simplify. We prefer to keep them in 
their general form but we will discuss later the numeri- 
cal differences between the complete and the systematic 
chiral expressions. The differences are one indication for 
the importance of higher-order chiral corrections. 

Although ImA is itself first order in isospin breaking 
we now make the usual (but scheme dependent) sepa- 
ration of the electroweak penguin contribution to ImA 
from the isospin breaking effects generated by other four- 
quark operators: 
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In order to perform such a separation within the CHPT 
approach, we need to identify the electroweak penguin 
contribution to a given low-energy coupling. In other 
words, we need a matching procedure between CHPT 
and the underlying theory of electroweak and strong in- 
teractions. Such a matching procedure is given here 
by working at leading order in 1/N C . Then, the elec- 
troweak LECs of 0(G 8 e 2 p n ) (n = 0,2) in ImA 2 lon ~ orap 
must be calculated by setting to zero the Wilson coeffi- 
cients C7, 6*8, Cg, C10 of electroweak penguin operators. 

Splitting off the electromagnetic penguin contribution 
to ImA in this way, we can now write e' in a more fa- 
miliar form as 
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The quantity fi e ff includes all effects to leading order in 
isospin breaking and it generalizes the more traditional 



parameter . Although fiie is in principle enhanced by 
the large ratio ReAp - 1 /Re A 2 the actual numerical anal- 
ysis shows all three terms in <|18H to be relevant when 
both strong and electromagnetic isospin violation are in- 
cluded. 

We present numerical results for the following two 
cases: 

i. We calculate £l c s and its components for a = 
(purely strong isospin violation). In this case, there 
is a clean separation of isospin violating effects 
in ImA2. We compare the lowest-order result of 
0(m M — vtid) with the full result of 0[(m u — md)p 2 ]- 

ii. Here we include electromagnetic corrections ex- 
plicitly, comparing again 0(m u — n%d, e 2 p°) with 
0[(m u — md)p 2 ,e 2 p 2 ]. In this case, the splitting 
between ImA™ 13 and lmA 2 on ~ crap is performed at 
leading order in 1/N C . 

In Table [I] the uncertainties in the "LO" entries are 
dominated by error propagation from the input parame- 
ters (these ranges account however for known 1/A C cor- 
rections). Apart from loops, new effects in "NLO" en- 
tries depend on ratios of next-to-leading to leading-order 
LECs. In Table |U we use the leading l/N c estimates for 
the ratios (GsN^/Gs, . . . , and estimate the uncertainty 
as discussed above. The final error for each of the quan- 
tities fiiB, Ao, /5/2, and fl e ff is obtained by adding in 
quadrature the LO error and the one associated to weak 



LECs at NLO. Moreover, only f 5/2 and ReA^/ReA 
depend on the ratio Gs/G 2 r. In these cases we rely on 
the phenomenological value implied by our fit 3]. Some 
of the errors in Table |U are manifestly correlated, e.g., in 
the LO column for a^0. 

The NLO results are obtained with the full expres- 
sions (|T3"j) , . . . , jlTST) . Using instead the simplified expres- 
sions corresponding to a fixed chiral order, the results are 
found to be well within the quoted error bars. The same 
is true for the alternative procedure of applying large- N c 
directly to the LECs GgNi, .... The isospin violating 
ratios in Table [I] are hardly changed at all in this case. 



1 to) 



4 





a = 
LO NLO 


Q / 

LO NLO 




11.7 15.9 ±4.5 


18.0 ± 6.5 22.7 ±7.6 


A 


-0.004 -0.41 ± 0.05 


8.7 ±3.0 8.3 ±3.6 


/s/2 





8.3 ±2.4 




11.7 16.3 ±4.5 


9.3 ±5.8 6.0 ±7.7 



TABLE I: Isospin violating corrections for e' in units of 10 -2 . 
The first two columns refer to strong isospin violation only 
(m u 7^ m,d), the last two contain the complete results includ- 
ing electromagnetic corrections. LO and NLO denote leading 
and next-to-leading orders in CHPT. 



Altogether, we expect our errors to be realistic estimates 
of higher-order effects in the chiral expansion. 

As first observed by Gardner and Valencia Q, there 
is a source of potentially large isospin violation through 
terms of the type (m u — m^GgiY, which get contribu- 
tions from strong LECs of 0(p 6 ) at leading order in large 
N c . The crucial LEC of <D(p 6 ) can be related to the 
mass splitting in the lightest scalar nonet that survives 
in the large- N c limit. The results in Tableware based on 
the evidence discussed in Ref. ^3 that the isotriplet and 
isodoublet states in the lightest scalar nonet are do (1450) 
and ^(1430), respectively (Scenario A). In the large- N c 
limit, the results are independent of the assignment of 
isosinglet scalar resonances. If instead the ao(980) is the 
isotriplet member of the lightest nonet surviving in the 
large-V c limit the couplings (m u — rridjGsNi would be 
strongly enhanced [J^, implying large (of the order of 
100 %) and negative values of f^iB- In this case, c'/e 
would be more than twice as big as in the isospin limit. 
We consider this an additional argument against such 
a scenario for the lightest scalar nonet (Scenario B in 
Ref. [13). 

Finally, we have investigated the impact of some sub- 
leading effects along the lines of Ref. |l2j. Although by 
no means a systematic expansion in 1/N C , those nonet 
breaking terms may furnish yet another indication for 
the intrinsic uncertainties of some of the LECs. The size 
of those terms depends on the assignment of isosinglet 
scalar resonances. Since nonet breaking effects are large 
in the scalar sector they affect most of the entries in Ta- 
ble [I] in a non- negligible way, although always within the 
quoted uncertainties. Employing again scenario A for 
the lightest scalar nonet [12j, i7 e ff in H17fl decreases from 
6.0- 10~ 2 to -1.4- 10~ 2 . 

We have performed a complete analysis of isospin viola- 
tion in e', due to both the light quark mass difference and 



electromagnetism. In particular, we have included for 
the first time isospin violation in the ratio ImAo/HeAo, 
parametrized by the quantity Ao. This ratio gets only 
a small contribution from strong isospin violation but 
the electromagnetic part is important, being dominated 
by electromagnetic penguin contributions. Both Ao and 
the purely electromagnetic AI = 5/2 contribution / 5 / 2 
interfere destructively with Qib to yield a final value 
^eff = (6.0 ± 7.7) • 10~ 2 for the overall measure of isospin 
violation in e'. If electromagnetic penguin contributions 
are included in theoretical calculations of ImAo/ReAo, 
Ao can be dropped in fl c g to a very good approximation. 
Finally, if all electromagnetic corrections are included in 
lmA /ReA , ImA 2 /ReA 2 , and Re^4 2 /ReA , fi e ff is es- 
sentially determined by f2iB- 
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